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Abstract 

We present preliminary measurements of CP- violating asymmetries for the decay B° — ► uK®. The 
data sample corresponds to 347 million BB pairs produced by e + e~ annihilation at the T(4S) 
resonance. For the decay B° — > luK®, we measure the time-dependent CP-violation parameters 
S = 0.62±g;§^ ± 0.02, and C = -0.43^;^ ± 0.03, where the first uncertainty is statistical and the 
second systematic. 
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1 INTRODUCTION 



Measurements of time-dependent CP asymmetries in B° meson decays through a Cabibbo-Kobaya- 
shi-Maskawa (CKM) favored b — » ccs amplitude ^ |2] have firmly established that CP symmetry 
is not conserved in the neutral B meson system. The effect, arising from the interference between 
mixing and decay proportional to the CP-violating phase (3 = arg (—V cc [V* b /VtdV t * b ) of the CKM 
mixing matrix manifests itself as an asymmetry in the time evolution of the B°B° pair. 

In the Standard Model, decays of B° mesons to charmless hadronic final states such as luK° 
proceed mostly via a single loop (penguin) amplitude with the same weak phase as the b — ► ccs 
transition M, but CKM-suppressed amplitudes and multiple particles in the loop introduce addi- 
tional weak phases whose contribution may not be negligible; see Refs. [51 El for early quantitative 
work in addressing the size of these effects. We define AS as the difference between the magnitude 
of the time-dependent CP-violating parameter S (given in detail below) measured in these decays 
and S = sin2/3 measured in decays to charmonium and a neutral kaon. For the decay B° — ► ujK°, 
these additional contributions are expected to give AS ~ 0.1 0|H], although this increase may 
be nullified when final-state interactions are included |Sj. A value of AS inconsistent with this 
expectation could be an indication of new physics [§]. 

We present an improved preliminary measurement of the time-dependent CP-violating asym- 
metry in the decay B° — > ljK®, previously reported by the BABAR and Belle Collaborations |1UI 111] • 
Charge-conjugate decay modes are implied throughout. 

2 THE BABAR DETECTOR AND DATASET 

The data were collected with the BABAR detector 12 at the PEP-II asymmetric-energy e + e~ 
collider. An integrated luminosity of 316 fb" 1 , corresponding to 347 million BB pairs, was recorded 
at the T(4S) resonance (center-of-mass energy y/s = 10.58 GeV). Charged particles are detected 
and their momenta measured by the combination of a silicon vertex tracker (SVT), consisting of 
five layers of double-sided detectors, and a 40-layer central drift chamber, both operating in a 
1.5 T axial magnetic field. Charged-particle identification is provided by the energy loss in the 
tracking devices and by the measured Cherenkov angle from an internally reflecting ring-imaging 
Cherenkov detector covering the central region. Photons and electrons are detected by a CsI(Tl) 
electromagnetic calorimeter. The instrumented flux return of the magnet allows discrimination of 
muons from pions. 

3 ANALYSIS METHOD 

From a B°B° pair produced in an T(AS) decay, we reconstruct one of the B mesons in the final 
state / = LoK®, a CP eigenstate with eigenvalue —1. For the time evolution measurement, we 
also identify (tag) the flavor (B° or B°) and reconstruct the decay vertex of the other B. The 
asymmetric beam configuration in the laboratory frame provides a boost of /?7 = 0.56 to the 
center-of-mass in the lab frame, which allows the determination of the proper decay time difference 
At = tf—ttag from the vertex separation of the two B meson candidates. Ignoring the At resolution 
(about 0.5 ps), the distribution of At is 

e -|At|/T 

F(At) = — [1 =f Aw ± (1 — 2w) (Ssm(Am d At) - Ccos(Am d At))]. (1) 
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The upper (lower) sign denotes a decay accompanied by a B° (B°) tag, r is the mean B° lifetime, 
Arrid is the mixing frequency, and the mistag parameters w and Aw are the average and differ- 
ence, respectively, of the probabilities that a true B° (B°) meson is mistagged as a B° (B°). The 
parameter C measures direct CP violation. 

The flavor-tagging algorithm 1 has seven mutually exclusive tagging categories of differing 
purities, including one for untagged events that we retain for yield determinations. The measured 
analyzing power, defined as efficiency times (1 — 2u>) 2 summed over all categories, is (30.4 ± 0.3)%, 
as determined from a large sample of B decays to fully reconstructed flavor eigenstates (-Bfl av )- 

We reconstruct a B meson candidate by combining a Kg with anw-> ir + ir~ir° candidate. We 
select Kg — > ir + ir~ decays by requiring the ir + 7r~ invariant mass to be within 12 MeV (~4<r) of the 
nominal K° mass and by requiring a flight length greater than three times its error. We require 
the 7r + 7r~7r° invariant mass (m^) to be between 735 and 825 MeV. Distributions from the data 
and from Monte Carlo (MC) simulations |T3] guide the choice of these selection criteria. We retain 
regions adequate to characterize the background as well as the signal for those quantities taken 
subsequently as observables for fitting. We also use the angle 9h, defined in the to rest frame as the 
angle of the direction of the boost from the B rest frame with respect to the normal to the uj decay 
plane. The quantity Ti = |cos##| is approximately flat for background decays and distributed as 
cos 2 Oh for signal decays. 

A B meson candidate is characterized kinematically by the energy-substituted mass toes = 



y (^s + p • pb) 2 /Eq — p 2 B and the energy difference AE = E* B — y/s/2, where (2?o,Po) and 

(Eb,Pb) are four- momenta of the T(AS) and the B candidate, respectively, and the asterisk denotes 
the center-of-mass rest frame. We require \AE\ < 0.2 GeV, 5.25 < toes — 5-29 GeV, |A£| < 20 ps 
and o"a< < 2.5 ps. 



use the angle 9t between the thrust axis of the B candidate and that of the rest of the tracks 
and neutral clusters in the event, calculated in the T(45') rest frame. The distribution of cos 9t is 
sharply peaked near ±1 for jet-like qq pairs and is nearly uniform for the isotropic B decays; we 
require |cos^t| < 0.9 . 

From MC simulations of B°B° and B + B~ events, we find evidence for a small (0.3% of the 
total sample) BB background contribution. We have therefore added a BB component to the fit 
described below. 

We use an unbinned, multivariate maximum-likelihood fit to extract signal yields and CP- 
violation parameters. We use the discriminating variables toeSj AE, m^, TL, and a Fisher dis- 
criminant T 14 j. The Fisher discriminant combines five variables: the polar angles with respect 
to the beam axis in the T(4S) frame of the B candidate momentum and of the B thrust axis; the 
tagging category; and the zeroth and second angular moments of the energy flow, excluding the B 
candidate, about the B thrust axis |14j . We use At to extract the CP- violation parameters, S and 



We define the probability density function (PDF) for each event i, hypothesis j (signal, BB 
background and qq background), and tagging category c: 




To help reject the dominant background from continuum e + e 



qq events (q 



n, d, s, c), we 



C. 



(2) 



where a\ t is the error on At for event i. We write the extended likelihood function as 



(3) 
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where Yj is the fitted yield of events of species j, fj >c is the fraction of events of species j for each 
category c, and N c is the number of events of category c in the sample. We fix / s i glC and f BB to 
/s flav ,c> the values measured with the large -Bfl av sample pQ. 

The PDF 7 3 sig(At, o"At,c) is given by F(At) (Eq. P) with tag category (c) dependent mistag 
parameters convolved with the signal resolution function (a sum of three Gaussians) determined 
from the -E>fl av sample. The other PDF forms are: the sum of two Gaussians for all signal shapes 
except Ti., and for the peaking component of the 77137,- background; the sum of three Gaussians for 
Vqq(At,c) and V BB (At,c); an asymmetric Gaussian with different widths below and above the 
peak for Vj{!F) (a small "tail" Gaussian is added for VqqiT)); Chebyshev functions of second to 
fourth order for the TL distribution for signal and the slowly- varying shapes of the AE, m^, and 7i 
distributions for backgrounds; and, for V q q(mEs), a phase-space- motivated empirical function [T5] . 
with a small Gaussian added for V BB {m^). 

We determine the PDF parameters from simulation for the signal and B B background com- 
ponents. We study large control samples of B — * Dn decays of similar topology to verify the 
simulated resolutions in AE and mEs> adjusting the PDFs to account for any differences found. 
For the qq background we use (rriEs , AE ) sideband data to obtain initial PDF-parameter values, 
but ultimately leave many of them free to vary in the final fit. 

4 RESULTS 

The free parameters in the fit are the following: the signal, BB background, and qq background 
yields; the three shape parameters of VqqiF); the slopes of V q q{AE) and Vqq{mj,^)\ the fraction 
of the peaking component of Vqqim^); the 7Bes background shape parameter £ |TS]; S; C; the 
fraction of background events in each tagging category; and the six primary parameters describing 
the At background shape. The parameters r and Amj are fixed to world-average values |16j . 
Tabled shows the results of the fit. The errors have been scaled by ~1.10 to account for a slight 
underestimate of the fit errors predicted by our simulations when the number of signal events is 
small. 

Table 1: Total sample size, detection efficiency, signal yield, BB background yield and CP- 
asymmetry parameters S and C from the fit. 



Quantity w -^s 

Total fit sample 
Eff. (%) 
Fit signal yield 
BB yield 
S 

c 



Fig. H shows projections onto the fit variables for a subset of the data (including 45-65% of 
signal events) for which the signal likelihood (computed without the variable plotted) exceeds a 
threshold that optimizes the sensitivity. Fig.[2]shows the At projections and asymmetry of the time- 
dependent fit applying the same event selection criteria as for Fig. ^ Based on explicit variation 
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Figure 1: B candidate projections for B° — ► uK® of (a) w-es, (b) A.E, (c) T, (d) "H, and (e) m^, 
shown for a signal-enhanced subset of the data (points with error bars), with the fit function (solid 
line), and the background components (dashed line) overlaid. 



of C with S allowed to float, we find the correlation between S and C to be negligible. 



11 



BaBar 
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Figure 2: Projections onto At for B — > W-ftTg, where ic"p is the decay time for the signal B meson. 
Data (points with errors), the fit function (solid line), background component (dashed line), and 
signal component (dotted line) are shown for events in which the tag meson is (a) B° and (b) B°, 
and the asymmetry (Ngo — N^ )/(N^o +N^ ) is shown in (c), where N indicates the total number 
of events passing the same cuts as for Fig. ^ 
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5 SYSTEMATIC UNCERTAINTIES 



We estimate systematic uncertainties in S and C from the following sources: potential dilution due 
to BB background (0.01); variation of the PDF shapes used in the fit (0.01); knowledge of the 
parameters used to model the signal At distribution (0.02); and interference between the CKM- 
suppressed b — > ucd amplitude and the favored b — > cud amplitude for some tag-side B decays 
|17j (0.02 for C, negligible for S), where the value in parentheses is the size of the estimated 
systematic uncertainty. By applying distortions to MC samples and refitting all tracks, we find 
that the uncertainty due to possible SVT misalignment and position and size of the beam spot 
are negligible. The uncertainties in the parameters of fits to the -Bg av sample are used for the 
uncertainties in the signal PDF parameters: At resolutions, tagging efficiencies, and mistag rates. 
Published measurements ^5] are used for tb and Am^. Summing all systematic uncertainties in 
quadrature, we obtain 0.02 for S and 0.03 for C. 

6 SUMMARY 

In conclusion, we have presented preliminary results for the time-dependent asymmetry parameters 
for the decay B° -» ojK% S = 0.62±g;|g ± 0.02 and C = -0.43±g;| ± 0.03, where the first 
uncertainty is statistical and the second systematic. If we fix C = 0, we find S = 0.63 j^jj, where 
the uncertainty is statistical only. This value of S and the world-average value of sin2/3 IZj yield 
a value of AS = —0.09 ± 0.31, in good agreement with the Standard Model expectation near zero. 
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